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Abstract 

We construct a measure- valued equivalent to the spatial A-Fleming-Viot process (SLFV) intro- 
duced in [Ethos ] . In contrast with the construction carried out in [EthOS] , we fix the realization of 
the sequence of reproduction events and obtain a quenched evolution of the local genetic diversities. 
To this end, we use a particle representation which highlights the role of the genealogies in the at- 
tribution of types (or alleles) to the individuals of the population. This construction also enables us 
to clarify the state-space of the SLFV and to derive several path properties of the measure- valued 
process as well as of the labeled trees describing the genealogical relations between a sample of 
individuals. We complement it with a look-down construction which provides a particle system 
whose empirical distribution at time i, seen as a process in i, has the law of the quenched SLFV. In 
all these results, the facts that we work with a fixed configuration of events and that reproduction 
occurs only locally in space introduce serious technical issues that are overcome by controlling the 
number of events occurring and of particles present in a given area over macroscopic time intervals. 
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Key words and phrases: random environment, generalised Fleming- Viot process, spatial coales- 
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1 Introduction and motivation 

The spatial A-Fleming-Viot process was introduced in |Eth08| and |BKE10| as a new framework for 
studying the genetic evolution of a population distributed over a continuous space. In these papers, 
two-dimensional spaces such as or the 2-dimensional torus were chosen since this seems natural for 
modeling the habitats of biological populations. The model, however, allows to consider any domain 
T C M*^ as a geographical space, and any compact metric space K of types (or 'alleles'). We refer to 
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|BEH09) . IBEVIO) . [BEVllj . |EV12) and |BEV12 ; for first studies of tliis process and its relations to 
other models. 

One noteworthy difference with previous models of structured populations is the following: while 
in the latter the reproduction mechanism can be phrased in terms of clocks carried by the individuals 
and deciding of reproduction and death times, in the spatial A-Fleming-Viot process it is driven 
by a Poisson point process of extinction/recolonization events (of various ranges r > and impacts 
u S (0, 1]). For brevity, we shall address such a configuration of events as an environment. Once an 
environment lo = {{ti, Zi,ri,Ui), z € /} has been fixed, local populations change only at the times ti 
of an event, and only within the corresponding ball M{zi,ri) C T. During such an event, a parent is 
chosen uniformly at random within ]B(zj,rj) and a fraction Ui of the local population there is replaced 
by offspring of that individual in such a way that the total mass of the (continuum) population in this 
region remains constant. These offspring inherit the allele of their parent, hence results a jump in the 
local allele frequency at that time. As we shall see, the Poisson point process formulation also enables 
us to model the ancestral process of a sample of individuals, which, thanks to an inherent assumption 
of neutrality, is a process independent of types: backwards-in-time, all ancestral lineages which are in 
the geographical area of an event are chosen to be (or not to be) affected by the event in an i.i.d. 
manner, that is by tossing a coin with success probability u. The affected lineages then coalesce into 
a single ancestor whose position is uniformly distributed over the area of the event, and the ancestral 
process starts again from this new configuration. We shall assume that the random environment is a 
Poisson process with a possibly infinite intensity. 

In a non-spatial situation (i.e. with T consisting of a single element only), the random environment 
reduces to a Poisson process of so-called u-mergers, with intensity measure dt®i'{du) = dt u~'^ A{du) . 
In |BLG03| . the resulting forward evolution of type frequencies was analysed as a measure- valued 
process and called A- Fleming- Viot process. 

The spatial A-Fleming-Viot process (SLFV, in short) was formally constructed in [B_EyiO] using 
a method of Evans [Eya97j. This more analytic approach, which we re-describe in Section [2.11 relies 
on a characterization of the semigroup of the SLFV through a family of duality relations with the 
genealogical process mentioned above. 

Our first motivation for this work is to relate the forwards and backwards evolutions in a more 
detailed manner, by first defining a reproduction model conditioned on the environment and then by 
attributing the types thanks to an explicit use of the genealogical relations between the individuals of 
the population. This construction, which leads to Theorem 1, provides a more probabilistic approach 
to the definition of the process and enables us to include a mutation mechanism (which was not 
incorporated in the analytic construction and in fact seems less tractable there). In order to further 
disentangle the different sources of randomness acting on our population, in Section 12.3.11 we extend 
the environment lo by adding a fifth coordinate recording the spatial location of the parent during 
each event. This allows to define a parental skeleton, out of which emerge all the genealogical trees 
relating a sample of individuals. In Section I2.3.2( we randomly label the individuals sitting at each 
site by numbers from [0, 1], and extend the environment by a sixth coordinate that records the label 
of the parent during each reproduction event (in the spirit of Bertoin and Le Gall's flows of bridges, 
see [BLG03]). We show that in the absence of mutation, conditionally on this extended environment 
the SLFV can be seen as a deterministic flow on the field of local allele frequencies which at any time 
conserves the geographical density of the population (here assumed to be constant over space), and 
transports types through an accumulation of local reproduction events. Mutations may distort this 
flow by acting independently on each ancestral line between the reproduction events. 

A second motivation consists in clarifying the state-space of the SLFV and its topological properties. 
This is done in Lemma 12.11 Under the basic integrability assumption of [BEVIOJ (see condition ([T]) 
below) we obtain some general properties about the process: Lemma [2. 61 shows that the quenched (i.e., 
for a fixed environment) SLFV without mutation has cadlag paths of finite variation with probability 
1; Proposition 15. 2l tells us that, when the fraction of individuals replaced during an event is always less 
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than 1, the genealogical process of a countable sample of individuals never comes down from infinity 
(see Section [5] for a definition of this property). 

Thirdly, we propose a representation of the SLFV as the empirical distribution of the paths of an 
infinite collection of individuals. This is done in Section |4] through a look-down construction, which 
enables us to approximate the quenched spatial A-Fleming-Viot process by the empirical distribution of 
a large but finite set of individuals, uniformly over compact time intervals. The proof of our Theorem 2 
is inspired by that of Theorem 1.1 in |BB+09| , which considers non-spatial A-Fleming-Viot processes. 
With the additional spatial structure, the fact that reproduction events happen only locally brings in 
new technical subtleties. 

Since their introduction by Donnelly and Kurtz in )DK96) and |DK99| . look-down constructions 
have proved useful in different contexts. For example, in |BB+05| it is used to show that the genealogy 
of an a-stable branching process (with a E (1,2)) can be described as a randomly time-changed (3- 
coalescent. In f PW06| . it is at the core of the study of the evolving coalescent encoding the genealogy of 
the whole population in a Wright-Fisher model. Finally, in | iLZ12j it is used to show that a given class 
of A-Fleming-Viot processes ('coming down from infinity' sufficiently fast) with Brownian mutation 
have compact support at any positive time. Such a construction was established for the non-spatial 
A-Fleming-Viot processes in f BB+09| . In |DE+00| , a look-down construction for continuum-sites 
stepping-stone models was carried out, in which the genealogy consists in Feller processes which coalesce 
instantly upon meeting. Because the corresponding genealogies are fairly different, so are the details 
of their proof. Furthermore, they only show the convergence of the finite-dimensional distributions 
of their look-down process, whereas we prove its convergence in the space of cadlag paths. To our 
knowledge, there are no other constructions in the flavour of [DK 96j in which the spatial structure of 
the population influences their reproduction. Note however that Etheridge and Kurtz |EK12| propose 
a look-down construction of the SLFV in the spirit of [KRllj . and use it in particular to show that 
the flnite-density population model introduced in |BEH09| converges towards the annealed SLFV as 
the density of individuals tends to inflnity. 

In a more applied point of view, the different constructions carried out in Sections [3] and |4] allow 
us to interpret a 'real' population with reasonably large density as a Poissonization of the measures 
describing the state of the theoretical 'infinite' population. That is, at any time t we can think of 
our 'real' population as being the result of the evolution of the individuals with level less than some 
quantity N in the look-down construction (with not too small) between times and t, which form 
a locally finite Poissonian cloud of trajectories in T x K. Furthermore, the definition of the SLFV 
in a fixed environment may serve in reconstructing the main features of the history of a population. 
Indeed, the genetic diversity observed at a given locus within a sample of individuals gives us some 
information about quantities such as the local rate of coalescence, which depend essentially on the 
annealed evolution. We refer to |BEKV12) for more on inference issues related to the SLFV. However, 
if we now consider a second locus in the genome, or even a fraction of this genome (data which are 
now accessible through modern sequencing tools), the genetic diversities observed at all these sites are 
correlated first and foremost by the fact that they all flow through the same environment uj. Hence, 
any information on a; brought by the analysis of the diversity at one locus yields some constraints on 
the genetic diversity at other loci. A good understanding of the quenched SLFV is thus necessary to 
get a handle on the correlations between different parts of the genome and to devise statistical methods 
based on the available data to detect the presence of atypically big events acting as major catastrophes, 
or of natural selection acting on a given gene, or of any other kind of deviations from the hypothesis 
of neutral evolution made in this work. 

The paper is laid out as follows. In Section [21 we describe the "kernel- valued" version of the spatial 
A-Fleming-Viot process introduced in [Eth08] and state our main result, namely the existence and 
uniqueness of the quenched SLFV with mutation taking values in the space of Radon measures on 
T X IC. We also present the two extensions of the environment discussed above. Theorem 1 is proved in 
Section [3] by constructing the two-parameter semigroup of transition probabilities conditionally on the 



3 



configuration oj of events, for almost every u. While this construction focusses on fixed times s < t, in 
Section [3] we materialize the dynamics of the evolution by a look-down construction directly defined for 
all times. Though these two steps could in fact be summed up into a single one, we chose to keep them 
separated to put forward the role of the Poisson structure and of the exchangeability of the population 
in the construction of the SLFV, and because the addition of levels is somehow an elaboration on the 
ideas used in the first construction. In Section [5l we use the same kind of arguments as in the previous 
sections to address the question of coming down from infinity for the genealogical process. 

2 The spatial A-Fleming-Viot process and its sources of randomness 
2.1 The kernel- valued spatial A-Fleming-Viot process 

Let us start by describing the first construction of the SLFV, carried out in |Eth08| and |BEV10| . As 
the geographical space in which the population evolves we fix a domain T C M"^, and we write K for 
the compact space of all possible types. For a locally compact metric space (with E = T x'K ox E = 
as the generic examples), let JV[[E) (resp., Mi{E)) denote the set of all nonnegative Radon (resp., 
probability) measures on E, and let H be the quotient of the space of all Lebesgue-measurable maps 
(or probability kernels) p : T — > A1i(K) by the equivalence relation 



In other words, two maps p and p' are in the same equivalence class iff the measures dx p{x, dn) and 
dxp'{x,dK) on T X IC are equal. We shall give a mathematically equivalent description of H at the 
beginning of the next subsection. 

To specify the dynamics, let /x be a cr-finite measure on (0, oo) and {I'r, r > 0} be a collection 
of probability measures on [0, 1] such that the map r i— )• z/,. is measurable with respect to p. Let 11 
be a Poisson point process on R x T x (0, oo) x [0, 1] with intensity measure dt ® dz ® p{dr)vr-{du) . 
The random point configuration 11 acts as a random environment in which the population evolves: if 
(t, z, r, u) € n, then at time t a reproduction event occurs within the closed ball M{z^ r) C T, that is 

• a location y is sampled uniformly at random in ]B(z,r) and a type k is chosen according to 
Pt-{y, dn) (equivalently, a type is sampled according to the mean type distribution in B(z, r) just 
before the event), 

• for every x € M[z,r), pt{x,dn) := (1 — u)pt-{x,dn) + u5k- 

The value of p outside the ball remains unchanged, that is pt{x^ dti) = pt-{x, dn) for every x ^ M(z, r). 

The existence and uniqueness of a Markov process evolving according to this dynamics is proven 
in Section 4 of j BEVlO] under the condition 



(This is shown there only for T = in fact, but the proof is identical in all dimensions and for any 
domain in M*^.) The process is characterized by a family of duality relations, based on the following 
very simple idea. In the absence of mutation, the type of the parent is transmitted to all its offspring. 
Thus, to know the types of a few individuals sampled at time t, it suffices to go back to their ancestors 
at some reference time (say, 0) in the past, and to check which types they carried. Now, some of the 
individuals alive at time t may share a common ancestor at time 0, which introduces some correlations 
between their types. As a consequence, the natural dual process to consider is that tracing the ancestral 
relations between a sample of individuals, from time t on and back into the past. 

Let us thus imagine what the ancestral process of a sample of individuals should look like under 
the prescribed evolution. Here we forget about types, and for any € N, we represent the genealogical 
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relations of a sample of k individuals (labelled by 1, . . . , A;) by a process {Ah)h>o with values in the set 
of all partitions of {1, . . . ,k} whose blocks are marked by an element of T. In words, h = corresponds 
to the time at which our individuals are sampled and for any h > 0, each block of Ah contains the labels 
of all of these individuals who share a common ancestor h units of time back in the past. The mark of 
a block records the spatial location of the corresponding ancestor at that time. When convenient, we 
shall write 



See Figure E] (with h = t-s), where k = A, Nh = 2, Bl = {1, 2, 3}, Bl = {4}. 

As a start, suppose that we sample a single individual and trace back where its ancestors were. 
Since forwards-in-time the population evolves only at the epochs of events of H, so does the spatial 
location of an ancestral lineage. Going back by a sufficient amount of time, one will encounter the 
first event (t', z' , r' , u') in the past in which our individual was not only in the area ]B(z', r') hit at that 
time, but was also part of the fraction u' of the local population replaced by offspring of the elected 
parent. (Below we shall argue that condition ([1]) guarantees that the time it takes back to this first 
event in the past is a.s. strictly positive.) Hence, at the time h = t — t' of this event, the ancestral line 
of our individual jumps to the location of its parent, which is by construction uniformly distributed 
over M(z', r'). Using the same reasoning, we can then find an earlier event (t", z'\ r", u") during which 
this parent was born. The ancestral lineage at that time jumps to a location uniformly distributed 
over M{z" ,r") and stays there until the time of the event in its past during which it was born, and so 
on. Now, observe that the sequence of events experienced by the ancestral lineage backwards-in-time 
is governed by the law of 11. Consequently, the rate at which the lineage jumps can be computed 
directly in terms of /i and Vr- once the radius r and the impact u of the event have been chosen, the 
volume of centres such that the current location of a lineage belongs to the range of the event is 
Vol(B(^|j,r) n T) and the probability that the lineage belongs to the set of newborns is u. Hence, the 
time to wait before the next jump is exponentially distributed with parameter 



and at that time ^* jumps to a new location uniformly distributed over the area of the corresponding 
event. Under the condition stated in ([T]), this jump rate is bounded by a constant independent of the 
current location of the lineage. 

Now think of several individuals sampled at distinct locations. When a past event comes up that 
covers at least one of them, each of the lineages within its range has probability u of being an offspring 
of the parent chosen, independently of the others. Then, all the lineages which were born in this 
event trace back to the same parent and therefore merge at that time into a single ancestor (i.e., the 
corresponding blocks of Ah merge into a single block) whose location is uniformly distributed over the 
area of the event. The other lineages, outside the ball or inside but not within the pool of offspring, 
remain unaffected. Because nothing happens between the events in H that affect the genealogy, the 
description of ^ as a Markovian system of coalescing blocks whose marks evolve like a family of 
correlated jump processes is complete. 

There remains to find an appropriate set of test functions. To this end, let C{E) (resp., Cc{E)) 
stand for the set of all continuous (resp. continuous and compactly supported) functions on the space 
E. Let also (m, /) denote the integral of the function / against the measure m. If /c G N, F E Cc(T'^) 
and , . . . , E C (K) , let us define the function Gg by 
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G^{Ki,...,Kk) := ]^5j(/«j) 



(3) 
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and the function ; gfi, . . . , 5^) € C(H) by 



hip; F, 91,- ■ ■ , 




A slight modification of Lemma 4.1 in |BEV10| shows that the family of functions of the form 
Ik{- F, gi, . . . , gi:), with k, F and gi,...,gk as above, is dense in C(H). As a consequence, they 
constitute a relevant set of test functions. 

We can now state the duality relations upon which the characterization of the kernel- valued SLFV 
depends. In Theorem 4.2 of jBEV10| . the SLFV is defined as the unique H-valued Hunt process {pt)t>o 
satisfying: for every po G E, t > 0, k G N, F e Cc(T'') and gi,.. . ,gk G C{K), 



[lkipt;F,gi,..., 




(5) 

where E denotes the expectation w.r. to the process {pt)t>0: E is the expectation w.r. to the process 
A, randomized over the environment (we shall be more precise on this point in the next section), and 

Pfc(x) := {({l},xi), ... ,{{k},Xk)}. (6) 

The idea behind ([5]) is precisely that expounded above: the type of an individual living at time t is 
given by the type distribution poiCu ') the location of its ancestor at time 0, and if several individuals 
have a common ancestor at time (i.e., Card(i?j) > 1 for some l), they all share the same type 
with law po{(,f,-)- 

2.2 The SLFV as a measure- valued process in a random environment 

Let A4\ be the space of all nonnegative Radon measures on T x K whose 'spatial' marginal is equal to 
Lebesgue measure A on T. Using a well-known disintegration theorem (see e.g. [Kal02j . p. 561), it is 
not difficult to show that H is in one-to-one correspondence with this correspondence being 

m{dx , dn) = dx p{x , dn) . (7) 

(Here and below, we sometimes shortly write dx instead of \{dx).) We shall endow Aix with the 
topology %, of vague convergence and the associated Borel cj-field, and shall use M.\ as the state space 
for the measure- valued SLFV. Using the correspondence ([7]), the function 1^ defined in reads as 

Ik{m;F,gi,...,gk) := (m®^F(g)Gg>, 

where Gg was defined in ([3]). The following lemma guarantees that the set of functions of the form If. 
constitutes a wide enough family of tests functions, and reveals in particular that the convergence of 
a sequence {pn)n>i of elements of H (as defined in |BEV10| ) is equivalent to the vague convergence of 
the measures m„(da;, d/t) := dx Pn{x, dK) as n — )■ cxd. 

Lemma 2.1. a) The space {A4x,Tv) is compact. 

h) For A; G N, F G Cc(T'^), gi,...,gk G C(K), the function m 1— t- Ik{m; F, gi, . . . , g^) is 
continuous on M\. 

c) The linear span of the set of constant functions and of functions of the form Ik{- ', F; gi, . . . , g^), 
ken, F e Cc{T^) and 51, . . . , c/fc G C{K), is dense in C{Mx)- 
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Proof of Lemma 12.11 a) Since the 7^-closedness is immediate from the definition of Ad\^ it suffices 
to show that any given sequence (m„) in TWa has a subsequence that converges in the vague topology 
Tv To see this, note that T x K is locaUy compact and separable, hence by Theorem A2.3 in |Kal02| . 
A^(T X K) endowed with the vague topology is Polish (and in particular complete). Therefore, again 
by the quoted theorem, (A^AjT^) is compact. 

h) This follows from the fact that for any locally compact separable space E and /c G N, the mapping 
m — )• m®^ is continuous w.r.to the vague topologies on M.{E) and M.{E^). 

c) This space is an algebra that separate points in A4x and contains the constants. Since A4x is 
compact by a), the Stone- Weierstrass theorem gives us the result. □ 

In Section HI we shall need to put a metric on the space D_\4^ [0, oo) of all cadlag paths with values 
in A4\. We thus recall the following result (see, e.g.. Section 1 in |DK96| ). 

Lemma 2.2. There exists a sequence {fn)n>i of uniformly hounded functions in Cc(T x K) which 
separates points in M\. Furthermore, if {gn)n>i is such a sequence, then 

oo ^ 

d{m,m') := ^ \ {m,gn) - {m',gn)\, m,m' G Mx 
n=l ^ 

defines a metric for the topology of vague convergence on M\, while 

/•oo 

A{{mt),{m't)) := e~*d{mt,m[) dt 
Jo 

is a metric for the topology of locally uniform convergence on D_m^[0, oo). 

Let us fix again a measure n on (0, oo) and a collection {v^, r > 0} of probability measures on [0, 1] 
satisfying ([1]). Let 11 be a Poisson point process with intensity measure dt®dz(^ ^{dr)vr{du). We write 
P for the distribution of 11, and note that P assigns full measure to the set $7 of point configurations 
(jO = {ti, Zi,ri,Ui)i£i on M X T X (0, oo) x [0, 1] with the property that tj ^ t^i for i ^ i' and that for all 
s < t G M and each bounded subset B of T, 

rfui < oo. 

i:s<ti<t, ZiGB 

Recall the description of the ancestral process A given in Section [2.11 For P-a.a. uj £ Q and t G M, let 
■^pi*(x) stand for the law of the system of coalescing jump processes describing the genealogy of k of 
'individuals' sampled at time t at the locations x = (xi, . . . ,Xk), conditionally on the environment uj. 
That is, under P'^'* • 
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The evolution of A starts at time /i = in pA;(x) 
and uses only the events (ti, Zi,ri,Ui) E lv such 
that ti < t. 

Whenever one or more lineages belong to the 
range of an event, each of the lineages within 
B(zj,rj) takes part in this event with proba- 
bility Ui or remains unaffected with probability 



1 



Ui, 



independently of each other. All those 



lineages which are affected merge into a single 
lineage whose location is uniformly distributed 
over M{zi,ri). Then A remains constant equal to 
its new value At-ti until the next event of oj in 
the past which hits at least one of the lineages 
and for which at least one of these lineages takes 
part in the merging. 



Xi 



X2 



X3 X4 



t-s 



Figure 1 



Condition ([T|) guarantees that for any given i € M, for P-a.e. environment uj, with probability 1 
no lineage in A has an accumulation point of jumps in the time interval [0, cx)). Hence, for P-a.all u 
we can define P'^'* , (for all A: G N and A^'^-a.all x G T'^) on the space T> of coalescing cadlag paths 
with values in T. We shall write P* , for the ioint distribution on 0, x V defined by 



(8) 



For a £ T> and h > 0, we write ah for the restriction of a to the time interval [0, h] and note that at 
describes a forest of coalescing paths. Indeed, if o by time h has coalesced to paths, then ah consists 
of Uh trees a^, i = 1, . . . , n/j. (Again, recall Figure [TJ) 

Remark 2.3. Since the law P is invariant by translation in the time coordinate, for any bounded 
measurable function F : D M, the quantity 



is independent of t and will hence be written as 



Ep.(x)[i^M)] 



nxv 



F{a)F^^^^){duj,da). 



(See, e.g.. Corollary \2.4\ ) 



Another ingredient we shall need is a mutation mechanism. Let {ICt)t>o be a Feller process with 
values in K, defined on some probability space (T>, T , Q). For every k £ K and every genealogical tree 
a (rooted in a single individual and) having n leaves at some time /i > 0, let us write 



ffi,...,c/„ G C{K) 



(9) 



to characterize the distribution of the types at the leaves when the root has type k and types evolve 
along the branches of a according to the mutation process /C (we assume that this evolution occurs 
independently along distinct subtrees emanating from the same vertex). 

The A^A-valued SLFV with mutation can now be defined (in its quenched version) as follows. 



8 



Theorem 1. For¥-a.all u there exists a unique A4.\-valued time-inhomogeneous Hunt process (Mt)^^^ 



whose two-parameter semigroup is characterized as follows: For every s<t = s + /i€ 
ken,Fe C,(T'=) and 51, • • • G C{K), 



m G M\, 



F(Xl,...,XA;)E' 



UJ,t 



1=1 



dxi ■ ■ ■ dxk, (10) 



where pA;(x) was defined in ^ and m = dx p{x, •) as in 



Note that the right hand side of (jlOp is weh-defined. Indeed, it foUows from the above described 
construction of A that, conditionahy on N^i, the law of the ancestral locations . . . , ^^'^ is absolutely 
continuous w.r. to Lebesgue measure. 

In Section [3] we shall prove a slightly stronger version of Theorem 1. Using the homogeneity of the 
Poisson distribution P, we shall then be able to conclude the following annealed version of Theorem 1 
(recall the notation Ep^,(x) from Remark I2.3p : 

Corollary 2.4. There exists a unique Aix-valued Hunt process {Mt)t>o such that for every m E M.\, 
t>0,keN,FG Cc{T^) and 51, . . . S C(K), 



E. 



M^\F^Gg[ 
F(xi,...,3;fc)Ep^(x) 



Nt 



t=l 



p{^l,dKi) ■ ■ ■ p{^f',dHiNt) 



dxi ■ ■ ■ dxk 



(11) 



Remark 2.5. When there are no mutations (i.e., when fC is the constant process) , we recover formula 
and hence the existence and uniqueness result from \BEV10^ . 

Theorem 1 shows that the quenched spatial A- Fleming- Viot process is a strong Markov process 
with cadlag paths. In the absence of mutation, it has even stronger path properties. 

Lemma 2.6. For P-a. e. environment uj, the quenched SLFV without mutation has paths of finite 
variation V^-a.s. 



Lemma 12.61 is proved at the end of Section 14.2 



2.3 Extending the environment 

Before establishing the main results of this work, we gather here two constructions based on particular 
extensions of the environment. The first one will prove useful in the next sections, whereas the main 
interest of the second one is to relate the SLFV to Bertoin and Legall's flows of bridges. As discussed in 
the introduction, they also highlight the different layers of randomness which appear in the construction 
of the spatial A-Fleming-Viot process. 



2.3.1 Recording the locations of the parents 

Recall from the paragraph above ([S]) that the probability measures Pp^*(x) defined P-a.s, for each 
fixed t. As a preparation of the proof of Theorem 1, we shall specify a construction of a "good version" 
of P'^'* which will serve simultaneously for all t G M. 

To this purpose, we mark each event (tj, Zi,ri,Ui) with a parental location yi, uniformly chosen from 
M{zi,ri), and denote the realizations of the resulting Poisson point process hy ijj = {{ti, Zi,ri,Ui,yi) : 
i G /}. Its distribution will, by a slight abuse of notation, again be denoted by P. 
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For a given -0, we first construct a random graph that codes the ancestral relationships between 
the parental individuals living at times ti at locations yi. To this purpose, let H := {Ha/) be a family of 
independent, uniformly on [0, 1] distributed random variables, indexed by the pairs i ^ i' . For ti' < ti, 
let us define 



Gi 



Thanks to the integrability condition ([1]), for P-almost all ■0 and all i, the random configuration of 
time points {tj/ : ti' < ti,Gii' = 1} is locally finite a.s. Thus, for all such ip, there is an ?^-measurable 
event of full probability on which for all i G I there is a uniquely defined index Tr{i) G I (depending 
on ip and Ti) such that 

*7r(j) = sup{tj/ : ti> < ti,Giii = 1}. 

We now decree that the most recent event which af- 
fected a parental individual with index i was the event 
with index vr(z). In other words, the location of the 
ancestral lineage of the individual that is parental in 
the i-ih. event, when viewed backwards in time, sits at 
location yi for the time span ti — t^^i^ and then jumps 
to . From there, it evolves as the ancestral process 
of the parental individual with index 7r(z) (with which 
it has merged at time ti — t^^^^i-^ ) . Doing this for alH € / 
gives us the genealogical tree of all the parents chosen 
during an event, which we call . 

Figure 2 

Next, let us tie in the (locations of the) ancestral lineage of an individual sampled at time t at location 
X by deciding when it joins the parental graph . To this purpose, we use, for given and x, a family 
T-L^ of independent, uniformly on [0, 1] distributed random variables {Hf )i^j that are independent of 
the collection and set 

:= l{Hf >l-Mi}^{a:eB(zi,r,)}- 

Again thanks to the integrability condition ([1]), for all x and for '. 
uration of time points {tj : = 1} is locally finite a.s. Hence, 
an "H^-measurable event of full probability on which for all t G I 
7r(t, x) G I such that 

tn{t,x) = sup{tj :ti<t,Gi= 1}. 

The location process '* of the ancestral lineage of an individual sampled at time t at location x is 
then constructed as follows. Viewed backwards in time, the lineage remains located at x down to time 
t-!r{t,x) =• til s-iid at that time jumps to the parental location y^ in the event Cj. From that moment on, 
it has merged with the ancestral line of this parent and follows its evolution in Q"^ . 

This construction extends to finitely or countably many individuals Ij, j G J, where Ij is drawn 
from a prescribed location xj at time t, and J is a (finite or countably infinite) index set. For this 
we use independent uniformly on [0, 1] distributed random variables -ffp^ and insert H^''^ in place of 
H^^ to define the ancestral location process of Ij. For x = (xi, . . . ,Xk), we then obtain the desired 
version of Pp^*(x) averaging the distribution of {Ati, ■ ■ ■ ,Axf) over {yi)i^j. Let us also note that 
in this way we obtain, for a Poisson point process <I> on T with Lebesgue intensity measure, an a.s. 
construction of A^'*, the locations of the ancestral lineages of individuals sampled at the points of ^ 
at time t. This will be used at the beginning of the proof of Theorem 1'. 



ti 



t 



7r(i) 



P- almost all ip, the random config- 
for all X and all such ■0, we have 
[ there is a uniquely defined index 
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2.3.2 Recording the labels of the parents 



This subsection, though not required for the remainder of the paper, is intended to connect the SLFV 
setting to the representation of generalized Fleming- Viot processes in terms of bridges by Bertoin and 
Le Gall |BLG03| . To this purpose we define another extension of the environment, by assigning to 
each event {ti, Zi,ri,Ui,yi) a parental label li, independently and uniformly drawn from [0,1]. The 
realizations of the resulting Poisson point process will be denoted by x = {{ti, Zi,ri,Ui,yi,li) : i £ I}, 
and its distribution will again be denoted by P. 

The labels li will be used to encode differently the additional randomness obtained through the 
auxiliary random variables Hai and Hf in the construction described in the previous section. More 
precisely, following the "flow of bridges" idea of Bertoin and Le Gall |!BLG03| , we shall attribute a label 
to each individual in the sample and use the parental labels 1^ to trace back the desired ancestries. As 
in |BLG03| we define the elementary bridge associated with (u, I) as 

bi,u{'^) ■= (1 -u)w + ul[ii]{w), <w <l. 

The parental label of an individual with label v that is overlapped by an event with impact u and 
parental label / is defined to be {bi^u)~^{v)- In other words, under the inverse mapping (6/.u)^^ the 
interval [(l — u)l, {l — u)l + u] is mapped back to /, whereas the intervals [0, (1 — u)/) and ((1 — u)Z + u, 1] 
are stretched to [0, ^) and [l, 1], respectively. 

For each x G T and s < t, the bridge B^^ is designed below in such a way that {B^ (v) is the 
individual label of the ancestor (of the individual with label v sampled at location x at time t) that 
either lived at time s, or was parental in the most recent event before t which affected the ancestral line 
of the individual sampled at location x at time t. Each event affecting x between s and t contributes 
a jump to Bgf. and at the same time shrinks the "old bridge" (see Figure [3]). Since the parental labels 
are a.s. distinct and since the locations of the parents are a.s. different from x, when looking at the 
population at site x forwards in time, what we are trying to encode here is an immigration of jumps 
gradually replacing the continuous part of B^ -/- (which in turn represents the population at time t whose 
ancestors at time s were already living at x). 

More concretely, the way how two elementary bridges bu^^i^ and b^^^i^, with 6^1 ,h older than b^^^i^, 
affect X is not the usual composition bu2,i2 ° bui,ii{w) = &M2,i2(^Mi,h(^))' but is given by 

bu2,l2 O bui,h{w) ■■= {I - U2)bu^,l^{w) + U2l[l2,l]{w) 

= (1 - n2)(l - Ul){w) + (1 - U2)uil[i^^i]W + U2l[l2,l]{w), 0<w<l. 

That is, a jump of size U2 is inserted at the argument I2, and the old bridge is shrunk by the factor 
1 — U2 (see Figure [3]) . Doing this for all the affecting events (in the right order) gives the bridge 

Blt{w):=[ n E n i^-UJ)]u^^k^]{w), 0<w<l. (12) 

\j:s<tj<t J i:s<ti<t \j:ti<tj<t J 

The quantity po := 0(1" ^j) the fraction in the population that is left over at time t from the 

j:s<tj<.t 

population residing at location x at time s, and the quantity := | Yi ~ uj) \ Ui is the fraction 

\j:U<tj<t J 

in the population that is left over from the re-colonizers in the event €{. In particular, in the absence of 
mutation (or if we only consider the family structure of the population at a given site), we see that the 
type (or family size) distribution conditioned on the extended environment x is a deterministic function 
of its initial value and of x- Furthermore, for a random variable V that is uniformly distributed on 
[0,1], the probability that the pre-image {B''^^)~^{V) equals li is pi, and the probability that it is a 
continuity point of B^^ is po- 
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Figure 3: Flow representation of the population sitting at a given site after no, 1 and 2 events. After 
the first event, all the individuals whose levels belong to the interval [a, b] are offspring of the parent 
(a.s. chosen somewhere else in the area of the event) with label I. After the second event, this family 
has shrunk to the interval [a', b'], while the new parent with label I' has given birth to all the individuals 
in the interval [c',d']. 

Based on this last point, let us show that the ancestry of a sample of individuals can be reconstructed 
by a deterministic procedure involving only x s^nd the distribution of the level Vj of the individuals that 
we consider. Indeed, on the event {B^^)~^{V) = li for some i indexing an event, we set r(,x, V) := ti, 
and proceed back from {ti,yi) and from the individual level li, this time using those of the extended 
events in x that affect the location yj. 

The integrability condition implies that for all x and P-a.a. % 

Ui < oo. 

i:s<ti<t, \zi—x\<ri 

Therefore, for aU v G [0, 1], the path h ^ {B^_f^ ^)-^ (v) , < h < t, is a.s. of finite variation (though 
with jumps that generically will occur densely in time). Let us also note that for each h G [0, t], 
conditionally on h > t — T{x, V), the random variable f)~^{V) is uniformly distributed on [0, 1] 

when randomized over x- 

We can carry out the same construction for a sequence {xj, Vj) with independent Vi,V2, . . ., and the 
locations Xj not necessarily distinct. In this way we obtain a (x, Vi, V2, . . .)-measurable configuration 
of coalescing ancestral lines starting in {xj,Vj)j>i at time t, whose x-component, when randomized 
over y, has the distribution P* 

3 The SLFV conditioned on the environment 

In this section, we prove a slightly stronger version of Theorem 1, which will establish the latter a 
fortiori: for P-almost all u £ we characterize a two-parameter semigroup ^{m,dm') of transition 
probabilities on M\, where Qg ^{m, •) is the distribution of a random element dyts,t{^) hi A4\, which 
is the result of the transformation of m due to the mutations and the events in u occurring in the time 
interval [s,t]. 

Theorem 1'. For F-almost every a; G there exists a unique collection Q^^^{m,dm') of transition 
probabilities on Aix such that for every s<tGM, A;gN, Cc{T^), gi,...,gk G C{K), and 
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(13) 



F{xi, . . . ,Xfc)E' 



LJ,t 



N,,, 



/5(C/l,c^'^l)•••/5(4^t^'^A^J 



dxi ■ ■ ■ dxk, 



where here again m{dx,dK) = dxp{x,dK). Furthermore, ¥{duj)-a.s. we have 

(i) for all s <t, Qg^{m, •) is a continuous function of m ^ M\, 

(ii) for all s (zM., t ^ Q^o t > s, is a strongly continuous operator. 

(Hi) (Q^t) ^ two-parameter semigroup. More precisely, for all q < s < t, and for every F and 
gi, . . . ,gk as before we have 

Ql,{m,dm') I Q%t{m',dm"){m"^'',F®Gg) = [ Qlt{m,dm"){m"®^ ,F G^) . (14) 



Remark 3.1. Stated in words, relation (jl3p reads as follows: Let SUti_/j^((m) he a random measure 
with distribution Q^f-{m, •). Then a random k-sample drawn at locations xi, . . . ,Xk from 9Jl(_/i^t(m) 
arises in three steps: (1) Take, under P'^'*^ the random ancestry A started in pfc(x), (2) sample the 
types Ki, . . . , K^f^ drawn at locations from m, and (3) run the mutation process starting 

from the types ki, . . . , ktv^ along the genealogy A forwards in time, starting at time s = t — h, and 
up to time t. The types sampled at locations xi, . . . ,Xk at time t then result from the forest-indexed 
mutation process evaluated at time h. 



Before we prove Theorem 1', let us show that Coronary 12.41 is a straightforward consequence of it. 

Proof of Corollary 12.41 For any initial condition m G M.\, let us set Qt{m,-) = J ¥{dijj)QQ^{m,-). 
Using (lisp . dH) and Fubini's theorem, we obtain that for eyery /c G N, F G Cc(T'^) and gi, . . . ,gk G 

C{K), 

Qt{m,dm'){m"^^ , F ® G^) 



E 



/ F{x^,. 



,Xfc)E* 



Nt 



Nt 



L = l 



p{il,dKi)--- p{if\dKNt] 



dxi ■ ■ ■ dxk 



[ F(xi,...,x,)E^^(,) / flQf^ Udji^t) 



dxi ■ ■ ■ dxk 



(15) 



Hence ([TT]) holds and the existence of a suitable collection {Mt, t > 0} of random variables is proven. 
Next, Lemma [2.ir c) shows that the duality relations (|lip are sufficient to guarantee the uniqueness in 
law of each Mt. 

There remains only to show that {Qt, t > 0} is a semigroup. Indeed, once this has be established. 
Theorem 1' (z) and {ii) show that {Qt, t > 0} has the Feller property, and hence defines a Hunt 
process. Now for < s < t let us observe that, because of ([1] 



Qt{m,-) = / P((icj)Q^,(m,. 



¥{doj) Ql,{m,dm')Qlt{m',-). 



Since u i— > Qq ^ is measurable with respect to the events between times and s, a; i— t- Q^^ is measurable 
with respect to the events between times s and t, and P-a.s. no event happens precisely at time s. 
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then due to the independence property of the Poisson distribution P and by Fubini's theorem we can 
rewrite the r.h.s. as 

[ F{duj) [ Q'^^{m,dm') [ F{duj')Q'^\{m' , ■). (16) 

Since the environment distribution P (and together with it also the distribution of the genealogies ran- 
domized over the environment) is invariant under time-shift, we can write that J^¥{duj')Q'^t{m' , •) = 
F{duj')QQ f_g{m' , •). Hence ([TB]) equals Qsi^n, dm')Qt-s{fn' , •), which completes the proof of the 
semigroup property. □ 

Proof of Theorem 1'. Existence of Q^i{m). In a first step, we consider a Poisson point measure $ 
on T with intensity measure A. We write <I> = 'Ylij^j^xj and construct, for a given environment ip = 

{ti, Zi,ri,Ui,yi)i^l, the coalescing ancestral location processes {Atf)j^j as described in Section [2.3.11 
We run this processes for h = t — s units of time into the past and sample a type at each location 

of an ancestor according to the type distribution p{(,'f^,-)- Finally, we run the mutation process 
described in the paragraph around ([9]) along each of the ancestral trees we obtained, independently 
of each other. In other words, assuming the type was attributed to the t-th ancestor, we use that 
allele as a starting point of the evolution along the genealogical tree rooted at this ancestor, and the 
result of the mutation process at time h gives us the type of each of the points xj € $ such that j 
belong to the block B^. In this way, we obtain a point measure A/" on T x K, whose projection onto 
the geographical space forms a Poisson point process whilst the types assigned to the different points 
are correlated through the genealogy (and over the geography). 

In the second step, we perform the just described procedure for countably many i.i.d. copies 
$i,$2) . . . of resulting in a sequence A/i,A/2, ... of random point measures on T x IC. A simple 
but crucial observation is that the construction of the point measures J\fn is equivalent in law to the 
following construction: 

(i) Enrich the environment with the family Ti and construct the parental skeleton as in Sec- 
tion EXH 

(i) Allocate a type to each of the parents by running the mutation process along Q'^ , from time s 
on. This gives us a labelled skeleton S. 

(ii) For every x G <!>„, if the ancestral location process back from x coalesces with Q"^ (say, at 
time Tx back into the past), then start an independent mutation process from the type at this 
point of the skeleton and run it for the time span r^. If the ancestral process does not coalesce 
with Q"^ (and thus remains at x), then sample a type from the measure p{x,-). In both cases, 
allocate the resulting type to the 'individual' located at x in the point process. 

This description shows that, conditionally on S, the point measures A/i, ^2, . . . are independent and 
identically distributed. From this it is not difficult to conclude that := ^{Afi + - ■ ■+Mn) converges 

a.s. in the vague topology towards a (random) limit 9Jls,t as n ^ oo. Indeed, let us write Bci^ x IfC) 
for the set of all bounded functions with compact support on T x IC. By the law of large numbers 
(applied conditionally under ip and iS), for every / € -Bc(T x K), the sequence {{^J\f^, f))n^N converges 
P'^'*-a.s. and in to a random variable L{f). We thus put, for every / € Bc{T x K), 

{ms,tJ) ■■=L{f). (17) 

Up to restricting this definition to a countable set of functions / of the form 1^^ , where the compact 
sets {Aj, j € N} form a basis of T x K, we can conclude that P^'*-a.s., ^s,t is a (random) nonnegative 
Radon measure and ^A/"" converges a.s. in the vague topology towards Ws,t as n — )• cxd. On the 
other hand, again by the law of large numbers, for each F € Cc(T) and / := F 1^, we have that 
{^M^,f) —7- (A, F) a.s. as n ^ oo. Hence, with probability one Ws,t belongs to A4x. 
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Defining Qg^{m, .) as the law of ^s,ty it remains to prove (jl3p . For this we proceed by proving the 
following claim: 

For every G N, F G Cc(T'=), c/i, . . . , 9fc G C{K), we have 



lim E' 



UJ,t 



F{xi, ...,Xk)E 



Ul,t 

Pk (x) 



(18) 
dxi • • • dxk- 



Indeed, writing A/"" = ^{x^,K^) can compute 



E' 



n 



7iv,7fe 



UJ,t 



.n 



as ra ^ oo (the sum is set to be zero if M"^ has less than k points, which happens only if Vol(T) < oo 
and has a probability tending to as n — )• oo in this case). Writing := $1 + • • • + <^„, we have 



vr 



7i7^---7^7fe J=l 



rr 



E' 



UJ,t 



F(x,,,...,x,,)n5.(^7.) 

7l7^---7^7fc J=l 

r 

71 7^- ■■7^7* J=l 



But by construction, conditionally on the distinct values of X^-^ , ■ ■ ■ , X^^ we have, again with h := t — s, 



E' 



E' 



p{il,dKi) ■ ■ ■ pii^"" ,dKNj 



=: ^^(X). 

Thus, using the /c-th moment formula for the Poisson point measure = we arrive at 



E' 



^ -^(-^^71' • • • '^7j^s(^7i' • • • '^7fc) 

[ F{xi, Xk)'^g{xi, . . . ,Xk)dxi - ■ ■ dxk + o(- 
Jjk ^ \n 



+ o(- 

■ n 



which proves the claimed equality (jlSp . 

Denoting the (compact) support of the function F by Supp(-F) and using the fact that the number 
of points of <I>" within the support of F is a Poisson random variable with parameter nVol(Supp(-F)), 
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we can write 



E' 



UJ.t 



|((n-W")^',F®Gg>| 



< 



< 



k 



n 



LSupp(F), 



An II 1 CnfeVol(Supp(F))fe 



= C ||F|Unil9jlloo Vol(Supp(F))'= 

for a constant C independent of all other parameters. We can therefore use dominated convergence, 
together with the fact that ^A/"" converges vaguely to VJlg^t with P'^'*-probability 1, to conclude that 



lim E' 

n— >-oo 



'UJ.t 



Combining the above with (|18p yields (jl3p . and the proof of existence is complete. 

Uniqueness. By Lemma l2.ir c). the equalities (jl3p for all functions of the form 1^ are sufhcient 
to ensure that there is at most one distribution on Mx which satisfies them for any fixed w, m, s, t. 
Hence, uniqueness holds. 

Continuity with respect to m. This is a direct consequence of (I13p . Indeed, knowing that k lineages 
are sampled from the compact support of F, one can truncate the distribution of (^^, . . . , ^^*) uniformly 
in {xi, . . . ,Xk) G Supp(F) and turn the truncated density into a continuous function with compact 
support so that, up to an arbitrarily small error term, the expectation in the right-hand side of (jlSp 
is an integral w.r. to m®^^ of some continuous and compactly supported function. Assume now that 
a sequence (r?T,„)„>i of elements of Ad\ converges vaguely towards m (and so m®^ converges vaguely 
towards m®^ for every k). By the definition of vague convergence, for every (xi, . . . € Supp(-F) 
the expectation in the right-hand side of (jlSp written with m„ converges to the same expectation with 
m as the initial type distribution. Using dominated convergence (and working 'up to an arbitrarily 
small correction term') we can conclude that 



lim 



Q^^^(m„,dm')(m'®^F«)Gg> = / Q^^((m, dm')(m'®^ F ® G, 



oo. 



By Lemma [2.1( c). this guarantees that Qsti''^n, ■) converges weakly to Q^f-{m, •) as n 

Strong continuity at s. We follow the construction in Subsection 12.3.11 Here the good u^s are the 
ones for which the parental skeleton exists (with a locally finite jump intensity along all its lineages) 
and for which for A-a. all. x the set of jump points {ti : Gf = 1} is locally finite on the time axis. 

By Lemma [2.1( c). all we need to show is that for all co with the just described property we have 



lim 

e->0+ 



(19) 



for every m G Mx, A; G N, F G Cc(T'=) and 5i, . . . , c/fe G C{K). 

Here again, the key tool is the set of equations (jl3p . We can couple the P'^'^^^ just as we did in 
Subsection 12 . 3 . 1 1 to see that for Lebesgue-almost all x = (xi, . . . ,Xk) G Supp(F), 

ei!JJ^pf(x) j™Ps during {s, s + e]] = 0. 

Together with the Feller property of the mutation process, we can thus conclude that for Lebesgue-a.e. 
X G Supp(F) 



lim E 

£-i>0 



n 
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This result being independent of the measure m taken to attribute the types ki, . . . , k^, (|13p and 
dominated convergence give us that (jl9p is satisfied for every m G P-a.s. 

Flow property. Let us give a conceptual proof, based on Remark 13.11 (which is stated below Theorem 
1'). For brevity, we write OJtg ,,(m) for a random element in A4x that has distribution Q'^ ^{m, •). With 

that notation, we have to show that WV^^{Wlg g{m)) = DJV^^{m) for P-almost all uj. By Lemma [2.H for 
this it is enough to show that the type distribution of a random fc-sample drawn at locations xi, • • • , Xj^ 
coincide for both random measures. We start with analysing these for 571^^^(971^ ,,(m)). In order to 

obtain the random types at the locations (,j_g, ■ ■ ■ ,S,f-s" as required in step (2) of Remark 13. 11 we have 
to go further down in the ancestry, now starting from the partition {({l},Ct_s)) • • • ; ({-^i-s}) Ct^V)}- 
By the Markov property of A under Pp^*^^) (^^^ Section 12.3. ip , this amounts to running A for a total 

time t — q. Consequently, we have to sample at the locations ^f_q, ■ ■ ■ ■,S,flq'^ from the measure m, and 
from the resulting types run the mutation process forwards in time between times q and s in order to 
obtain the types at time s. These are used as the input in step (3) of Remark 13. 11 that is, as the initial 
conditions for another go of the mutation process between times s and t. By the Markov property 
of the mutation process /C, this amounts to running K, along the genealogy A between times q and t. 
In total, we have thus arrived at the types of a random fc-sample drawn at locations xi, . . . from 
97tg^^(m). Notice that this scheme being independent of the measure m, it works simultaneously for all 
m's. 

The proof of Theorem 1' is now complete. □ 

4 A look-down representation of the SLFV-process 

In the previous section we constructed the two-parameter semigroup (Qs t) given the configuration u 
of events. This was done, first for fixed times s < t, on top of the random genealogy of a sample whose 
locations had a Poisson(cA)-distribution (with c — )• oo). One ingredient for constructing this random 
genealogy was the process of parental locations in the events, which was independent of the sampling 
locations. In this section, we shall again work with a Poisson system of sampling locations "with infinite 
spatial density" but now this system will evolve in time, with the locations of the sampling and the 
parental locations in the events being coupled. That is, we shall always choose the parent among the 
individuals of our Poisson system. 

Let us start with an informal description. At time t = we start with a Poisson configuration 
of particles on T x [0, cc) with intensity measure \{dx) ® M. The first component is the particle's 
location, the second will be called the particle's level. While the levels stay fixed in time, the locations 
Cj(t), given oj = {ti, Zi,ri,Ui)i^i G fi, perform independent jump processes: at each time ti such that 
||Ci(^i~) ~ ^ ^ii the particle at level ij tosses a coin with success probability Ui, independently 
of everything else. If this coin comes up with "success", the particle jumps to a location Cji^i) that is 
chosen uniformly from M{zi,ri), again independently of everything else. 

Because of the integrability condition ([1]), the jump times of every single particle are a.s. locally 
finite, and because of the Poisson colouring theorem, the process (Ci(*))^j)j remains Poisson, with 
invariant intensity measure \(dx) ^ d£ (as we shall see below). Among all the particles that are 
affected by the event at time ti there will therefore be a.s. one whose level is smallest, say The 
key idea is now to declare the location Cj{i){ti~) as the parental location in the event at time ij, and 
to decree that all those (ti,Q{ti)) for which Q jumped at time ti are children of {ti — ,(j{i)iti~))- In 
this way, a genealogy is filled into the space-time point configuration (ti, Cj(^j)) in a look-down manner: 
certain ones of the particles at higher levels look down at a particle at lower level and copy its type. 

Let us now proceed with a more formal definition of the particle system in order to prove a look- 
down representation of the SLFV as stated in Theorem 2. We shall discuss possible generalizations of 
this construction in Section [4.31 
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4.1 A particle system and its look-down genealogy 

Let us define the forwards-in-time motion {Ct)t>o of a single individual in the environment a; S as 
follows: whenever lies within the range of an event (tj, Zj, r^, Uj) G u, the process at time ij does 
nothing with probability 1 — Ui, or jumps to a new location uniformly distributed over B(zj,rj) with 
probability Ui. We write P!^ for the probability measure on D under which C starts at x £ T (recall 
that P is the space of coalescing cadlag paths with values in T, defined just above ([8])). 

Let us observe that Lebesgue measure is reversible for the evolution of Indeed, writing Ki{x, dy) 
for the transition kernel of C, during the event of oj labelled by i, we have 

(1 - Ui)5^ + yQ^(]g(^.^^.^^ dy 

so that it is easy to check that 

dx Ki{x,dy) = dy Ki{y,dx), x,yGT. 

In particular, Lebesgue measure is conserved by the flow through the countably many events of lo 
(which is reminiscent of the fact - proved in the previous section - that 9Jt^((m), m £ Mx, has 
Lebesgue measure as a 'spatial' marginal for any s and t). 

Let us now define a Poisson point process A/" on T x P x [0, oo) with intensity measure dx P^{d(^)f^dL 
That is, we fix a Poisson point process {xj)j^j of locations at time 0, and launch a path from the 
point Xj. In addition, each path is labeled by a level in [0, oo) which we shall use to indicate who 
reproduces during an event. Note that, by the invariance of Lebesgue measure under the dynamics of 
the spatial distribution of the population is conserved at any time. That is: 

Lemma 4.1. For every t >0 and every c G [0,oo], the set {{Ci^^j) ■ i^jiC'^i^j) ^ -^i — c} forms 
a Poisson point process on T x [0, c] with intensity measure dx l{£<c}d^- 

Since the epochs of a reproduction event are deterministic when uj is fixed, we can in particular 
use Lemma [4.11 at the time tj of an event. 

The process M encodes which individual takes part in a given reproduction event: if the time of 
this event is > 0, then all individuals such that (^^__ ^ C,l_ are affected by the event, while the others 
remain unaffected. Let us now define the ancestral lineage of the individual (i, Ct ) by first tracing back 
C,^ to its most recent jump before t. If this jump happened at time ti-^, say, and if is the index of 
the path that had the lowest level among all those affected at time tj*, then the ancestral location 
remains at C,l for < h < t — ti^: and jumps to C^* at time h = t — ti^,. Back from ti^,^ we then proceed 
inductively. 

As we did for the ancestral process A (defined in Section [2]), for any finite set of distinct indices 
ji, ■ ■ ■ ,jk we call Vl := Vl{ji, . . . ,jk) the marked partition defined as 

={(4,Xi),...,(vr^,X^)}, 

where each vr^ contains the labels of all the individuals in {ji, ■ ■ ■ ,jk} at time t who share a common 
ancestor at time t — h, the second component stands for the spatial location of this ancestor at time 
t — h, and L/j denotes the number of distinct ancestors at that time. Notice that, since the Poisson 
point process on T x 2? x [0, oo) has infinite intensity measure dxV^ ^ di, it always contains infinitely 
many individuals even though T may have finite volume. 

The following lemma is crucial for the look-down representation of the SLFV. 

Lemma 4.2. For any finite set {ji, ■ ■ ■ , jk} , (J^h)he[o,t] ^c-^ same law as the ancestral process 
{A.h)he[o,t] under F^'* (recall the notation pk from 

PkiCt '■■■iCt ) 



. ^, (20) 
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Proof. By Lemma l4.ll at any time t >0 the set {(Ct)^j) • {^jiQ'^i^j) ^ ^} forms a Poisson point 
process on T x [0, oo) with intensity measure dx ® di. Furthermore, the environment u is fixed and so 
the times tj at which the events take place are deterministic. For these two reasons, we first claim that 

Claim 4.3. During any given event, the spatial location just before the event of the affected individual 
with lowest level is uniformly distributed over the range of the event. 

Indeed, let us write (tj, Zi,ri,Ui) for this event. By the thinning property of Poisson point processes, 
the set of affected individuals forms a Poisson point process on M(zi,ri) x [0, oo) with intensity measure 
Uil^(^Zi,ri)ix) dx (8> di, and thus levels and locations are attributed independently. The result is then 
straightforward. 

Second, let us define r as the quantity in [0, t] such that t — r is the most recent time at which one 
of our k paths jumps before time t; the time t — r is necessarily the epoch ti of a reproduction event. 
Let us show that r has the same distribution as the first time at which A jumps under PpJ*, where we 
wrote instead of PkiCt^ , ■ ■ ■ , Ct'') fo notation. For u = 1, . . . ,k, let ajy(tj) stand for the set 

of indices of the events of uj occurring in the time interval (t j , t] and whose range overlaps (^j" . Let us 
now observe that, by independence of the paths l^^ (which jump or remain unaffected independently 
of each other during an event), we can write 

F^[t — T < ti] = [none of the k paths is affected during (ti, t]] 

k 

= np"[cf = cr, ys€{t„t]] 

u=l 
k 

=n n (i-^'^)- (21) 

But coming back to the definition of A under P"^'*^ we see that (j2ip is exactly the probability that 
none of the lineages starting from locations Ql^ , . . . , C,l^ jumps before time t — ti. Hence, the first time 
r at which V jumps has the same law as the first time at which A jumps. 

Thirdly, again by the independence oi C,^'^ ,. . . , Q^'' , all individuals present in the area of an event 
just before (and therefore just after) the event occurs choose independently of each other whether they 
take part in the event and jump, or not. Therefore, during an event affecting at least one of them, the 
choice of who is affected or not is decided by a set of independent Bernoulli r.v.'s, as in the evolution 
of A. By Claim the location of the individual with lowest level, onto which all affected individuals 
look down to find an ancestor, is uniformly distributed over the area of the event. Furthermore, the 
location of an affected individual is resampled independently of its current position during the event, 
so that the parental location is independent of the position of its offspring. This is precisely what 
happens to A during a merger. 

Lastly, one can pursue the analysis by defining the most recent time t — T2 such that the paths 
(^■^1, . . . , d^^" do not jump during the time interval (t — T2,t — t). The same reasoning gives us that T2 
has the same distribution as the epoch of the second jump of A under Ppf , that the set of individuals 
who look down on the lowest level during the event occurring at time t — T2 has the same law as the 
second merger of A and that the location of the ancestor is uniformly distributed over the area of the 
event. Carrying on in this way and using the facts that both V and A are finite-rate jump processes, 
we can conclude that they are equal in distribution. □ 

4.2 Look-down representation of the SLFV-process 

As in Section [HI let us now consider a Feller process {ICt)t>o with values in the compact type space K. 
Let B stand for the generator of /C. For convenience, we shall see B as an operator on Cc(T x K) that 
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acts only on the second coordinate. With this in mind, from now on we assume that the domain T>[B) 
oiB is dense in Cc(T x K). 

Using the particle system together with its genealogy from the previous subsection, for any measure 
m € M\, we construct a configuration of T x K- valued paths ((^f,JCl,ij)j£j by choosing the /Cq 
independently with distribution p{xj, ■), j J, and letting the IC^ evolve according to the look-down 
genealogy of the particle system. That is, mutations occur independently on each level and whenever 
a path j looks down on a path i with lower level, IC^ jumps to the current value of /C* at that time. 

For any n G N and t > 0, let us define 

j:e,<n 

From Lemma 14.21 and the fact that at any time, the set of locations of the particles forms a Poisson 
point process with Lebesgue intensity, we obtain that for P-almost every uj 

under P"^, where the marked point configuration AA" was defined in the proof of Theorem 1'. By an 
argument involving exchangeability and thinning, this extends, for any fixed time t > and all n € N, 
to the distributional equality 

{M},2Ml nMn = {M\M\ . . . , AA"). 

Hence, for any t > 0, the sequence (M")„>o has the same distribution as the sequence ^A/"", which, 
as was shown in the proof of Theorem 1', converges P*^ a.s. for P-almost every uj. Thus, for every 
/ G V{B) n Cc(T X K), one can define the a.s. hmit 

{Mr J) := hm (Mr,/) = lim i f{QiM)- (22) 

j:l.j<n 

This convergence holds in fact in a pathwise manner, as shown by the following result. Recall from 
Lemma 12.21 the topology of uniform convergence over compact time intervals with which L'_a/(^ [0, oo) 
is equipped. 

Theorem 2. ForV-a.e. environment uj , the sequence (M")„>i converges V^-a.s. towards the process 
{Mr)t>o, uniformly over compact time intervals. Furthermore, M°° has the same law as the quenched 
spatial A- Fleming- Viot process of Theorem 1 with initial condition Mq = m. 

Our proof of Theorem 14. 21 will be guided by the proof of Theorem 1.1 in |BB+09| . The new difficulty 
which arises in our setting is due to the need to control the number and the spatial distribution of 
all the particles with level less than n lying within the compact support Sf (for the 'geographical' 
coordinate) of some function / G Cc(T x K) in a given (deterministic or random) time interval. This 
is achieved through the following lemma. 

Lemma 4.4. For every cj G and every T >0, let 6 = 9{uj,T) be defined by 

9:= [ P^((Ct)o<t<Tn5//0)dx. 

JT 

Then, is finite with ¥ -probability 1. 

Proof of Lemma 14.41 This property is obvious when T has finite Lebesgue measure, and so let us 
concentrate on the case where T is unbounded. By Fubini's theorem, we have 

E[0]= / E[P^((Ct)o<t<TnS//0)]dx= / P4(G)o<t<T n 5/ / 0] dx, 

JT JT 
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where (in analogy to the notation introduced in Remark I2.3p := J P^{.)¥{du}) denotes the annealed 
probability measure of C- But under P^;, <^ is a finite-rate jump process whose instantaneous jump 
rates are bounded uniformly in the location by the quantity in ([1]). Let us assume that the jump rate 
i7 is independent of the location and equal to this quantity (as, e.g., when T = W^). Indeed, boundary 
conditions only slow down the evolution of a lineage by making the balls M{x,r) C T smaller than the 
d-dimensional ball B(x, r) C W^, and so the desired property will remain satisfied even when the jump 
rate of C is inhomogeneous in space. To simplify the notation, let us also suppose that Sf = M{0, a) for 
some a > 0. Again, we do not loose generality with this assumption since there exists a > such that 
Sf C ]B(0, a), and the probability of entering Sf before time T is bounded by that of entering B(0,a). 
Since the volume of possible centres for an event of radius r overlapping both x G T and a subset of 
Sf is bounded by the volume Vr of the ball B(a;,r), we can write that 



1 /"°° / \ VolfjS ) 

Px[C liss in Sf after its first jump] < — / ^( / uUi~{du) ) — 



fi{dr) 



(Note that a ball overlapping both x and Sf has radius at least (|x| — a)/2, then the volume of possible 
centres is bounded by Vr, the lineage jumps with probability u chosen according to Vr{du) and finally 
the probability that it jumps into Sf is bounded by Yo\{S f) /Vr). For x G B(0,3a), we shall simply 
bound this probability by 1. Integrating against Lebesgue measure, turning to polar coordinates, and 
finally using Fubini's theorem, we obtain that 

/ P^: [C lies in 5/ after its first jump] dx 
Jt 

< Vol(B(0, 3a)) + ^"^^^Z)*^^ Rd-i uuridu)^ fiidr)dR 



2 

< Vol(B(0,3a)) + '^^ (^j\ur{du)^ j^"^^" R'^-^ dR fi{dr) 

/"OO rl 

< Vol(B(0, 3a)) + C / r'^ uvr{du) ii{dr) := C < oo 



a JO 



by ([T]), where is a constant depending only on the dimension d, and C, C > depend on d and /. 

As concerns the subsequent jumps, recall from the beginning of Section [4. II that Lebesgue measure 
is invariant under the evolution of As a consequence, calling r the first time at which jumps and 
using the strong Markov property of we have that 

/ Px[C liss in Sf after its 2nd jump] dx = E^[P^^ [(^ lies in Sf after its first jump]] dx 
Jt Jt 

= / Px [C liss in after its first jump] dx < C . 
Jt 

Proceeding by recursion, we obtain that the integral over T of the probability that C, belongs to Sf 
after its i-th jump is bounded by C for any i>l. But by our assumptions, the number of jumps made 
by in the time interval [0, T\ is independent of its starting point and is a Poisson r.v. with parameter 
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JT. Hence, since the number of jumps of Q is also independent of its trajectory we can write that 

no]= [ Px[{Ct)o<t<TnSfy^(D]dx 
Jt 

Px [k jumps in [0, T]] P^; belongs to Sf initially or after one of its first k jumps] dx 

) 

oo / - ^ k 



'^k=0 



< e~"^'^ — — — >^ / Pa; \C lies in Sr after its ith jump] dx 



k=0 i=0 

< g- W_l_ (Vol(5/) + kC) = Vol(5/) + CJT < oo. 

k=0 

We can thus conclude that 6 < oo for a.e. environment. □ 

We can now turn to the proof of Theorem 2. 

Proof of Theorem 2. Let us start by the first statement. The key ingredient of our proof is adapted 
from Lemma 3.2 in |BB+09| : 

Lemma 4.5. For P-a.e. environment oj, for every T,e > and every f G T^{B) PI Cc(T x K), there 
exists a summahle sequence ((5„)„>i such that for every n € N, 



sup |(M,",/) - (M^,/)! > lle\ < 6n. 

\0<t<T J 



Indeed, let us fix / G ^^(-B) fl Cc(T x K) and for every s > 0, let T<j be defined by 

\t<s}u\o\{^{z,r)r\Sf), (23) 

(t,2:,r,u)6ij 

where here again SfdY stands for the compact support of / (for the 'geographical' coordinate). This 
quantity will give us an upper bound on the sum of the jumps of (M", /) over the time interval [0, T\. 
Note that 

E[Tt] = / / / / uNo\i^{z,r)r\Sf)vr{du)dzii{dr)dt 

Jo Jo JSf+B{0,r) Jo 

<T / u{Cir'^ AVol{Sf))ur{du)dzfi{dr) 

Jo JSf+B(0,r) Jo 

<TC2 / ur'^ur{du)fi{dr) +TCi / uYol{S f +M{0,r))ur{du)fi{dr) 
Jo Jo Jcz Jo 

oo fl 

d. 



< TC5 / / ur'^i^r{du)^i{dr) < 00 (24) 
Jo Jo 

for some constants Ci, C2, C3, C4, C5 > 0, where the last line uses the integrability condition (JT]). As 
a consequence, with P-probability 1 the quantity is finite. 
Now, let us define the deterministic times (ai)i>i by 

ai := inf It: > ^ I A 



and 

Oj+i := inf <^ t : Tt > T^^ + — I A ( + — 
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This sequence of times decomposes the interval [0, T] into at most := 2(Tj' +T)n^ subintervals, over 
which we shall control the fluctuations of M"" and M°°. To this end, let us define a sequence (di)i>i 
of stopping times (with respect to the filtration {a{{M^, s < t}), t > 0}) by 

a, := inf {t > a, : \{M^ J) - (M-,/)| > 6e}. 

Let us set 

:= |(M-,/) - {MIJ)\ V KM?^,/) - (M? ,/)|. 

Recall the quantity 9 = 0{u},T) defined in Lemma |4.4( and that the number of paths with level less 
than n and passing through Sf within the time interval [0, T] is a Poisson r.v. with parameter n9. By 
Lemma l4.4| 9 < oo for P-a.e. environments. In addition, 9 > Vol(5/) > P-a.s., which guarantees 
that n9 is a.s. commensurate with n. 

Let us now observe that, since /((^, k) = for any ( ^ Sf, for every t S [0, T] we can write 

jeOn ^ "'^ ie0„ 

where 0„ is the set of all indices j such that £j < n and (Ct )o<t<T H 5/ 7^ 0. The analysis of the 
parameter n9 made just above shows that, if we define the event An by 



An := <^ Card(e„) i 



— ,2n9 
2 ' 



(25) 



then there exists Cq > such that for every n G N, P'^(A„) < e But since 9 > 0, Lemma 3.1 in 

|BB+09| and the fact that Card(e„)/n becomes concentrated around 9 yield 



P'^fmaxHi >e) < P'^(A„) + ^P^(i7i > e and A^J < e"'^''' + 16(Tt + r)n2e-^^" 



(26) 



for some constant C7 > depending only on / and e. This controls the distance between M*^ and 
M°° at some discrete times. Since the definition of on bounds the variations of (M°°, /) over the small 
time intervals of interest, there remains to show that (M",/) does not fiuctuate too much over these 
intervals. 

For any % G {!,... and any index j € N, let Tij be the first time after at which the type 
and location of the individual at level £j is updated during a reproduction event, with the result that 
the value of f{C,^ ,1C^) changes (so that an individual outside Sf may be updated several times before 
the time Tjj at which it enters Sf). Splitting the evolution of the individuals of M" into its two 
components, mutation and reproduction, we can write that for any t G [ai,ai+i] 



+ {/(c^^^') - /(Ci..,,)->^W,)-)}- (27) 



n 



The second term in the r.h.s. of (j27p is bounded by A/"^(ai, 0^+1), where A/'r((aj, a^+i) denotes 

the number of particles j living at time a^, with level at most n and such that Tjj- < Oj+i. These 
particles are of three kinds: 

(a) either they belong to Sf at time ai and are affected by an event (s, z, r, u) G uj before time aj+i 
(each with probability u during this event). 
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(6) or they do not belong to Sj at time but the first event {t, z, r, u) that affects them happens 
before time Oj+i and brings them within Sj (each with probability nVol(S'/nB(2;, r))/Vol(B(z, r)) 
during this event), 

(c) or else they do not belong to Sj at time Oi, they jump at least twice between ai and aj+i and 
one of this jump (except the first one) brings them into Sf. 

We thus need to bound the number of each kind of particles, which we denote respectively by , 
M^l and M^j. To simplify the notation, let us also set e' := e/(2||/||oo)- 

First, by Lemma l4.ll the collection {(C«i)^j)) j £ J} forms a Poisson point process on T x [0, oo) 
with intensity measure dx (8) di. Hence, if all the particles initially (i.e., at time Oi) present in Sf were 
frozen until an event {s,z,r,u), the number of those affected by this event and with level less than 
n would be a Poisson r.v. Vs with parameter nuVol{Sf riM(z,r)). But some of these particles may 
have already been affected by a previous event, so that in fact Vg gives a stochastic upper bound on 
the number of particles affected for the first time and by this event. Considering now all the events 
occurring between and aj+i, we obtain that 

AA^J ^ Ps= Poisson r nuYol{SfnM{z,r))\ 

= Poisson (n(Ta,+i - T^J), 

where the equality is distribution comes from the fact that the Vs's can be chosen independent since 
they correspond to thinnings of the set {(Cai)^j) • C^i ^ Sf, ij < n} that are independent of each 
others. Now, by construction Tq,._^^ — Tq- < and so there exists a constant Cs > such that 

P'^(AA5"^ > ne'/S) < Prob(Poisson(l/n) > ne'/S) < e"^*'". (28) 

Using the same reasoning, the number of particles originally outside Sf and affected for the first 
time during an event (s, z, r, u) that brings them within Sf is stochastically bounded by a Poisson r.v. 
with parameter 

nuVol(B(z,r) \ S f))^^^^(^^^^ < nuYol{SfnM{z,r)). 
Hence, here again we have that 

P"'{^fn} > ne'/3) < e~^«". (29) 
Finally, Afj^^j is stochastically bounded by a Poisson r.v. with parameter 

„ / P^{at least 2 jumps i„ [a„a. + „-] and parses through S, in this interval) 4.. (30) 

JT 

Observe that is a stopping time with respect to the Poisson point process of reproduction events. 
Hence, using the Markov inequality and then Fubini's theorem, we can write that 

/ PS (at leas. 2 jumps in [a..., + „-] and passes through 5, in this interval) > .'J 

Jt 

< — Px(at least 2 jumps in [ai,ai + n~^] and passes through Sf in this interval) dx. 

Proceeding as in the proof of Lemma [4. 4| we obtain that the above integral is bounded by a constant 
times n~^. Consequently, 



3i < in ■ / Px{ — 2 jumps in [ai,ai+n ^] and passes through Sf in this interval) dx > e'/6 



<% 
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(Recall that is proportional to n?.) By the Borel-Cantelli lemma, we can conclude that for P-a.e. 

environments, there exists rig (a;) such that for every n > no, 



/ P^( > 2 ,.mps m [o„o. +„-] and parses through S, m this inte™l) < e'/O. 

Jt 



Vi < ir 



Assuming that our fixed environment belongs to this set of full P- probability and recalling (|3Up . we 
obtain that for every n > no(w) 



F^{3i < i„ : Mi'} > ne'/3) < Cion^e"^"". 
Combining ([28]), ([29]), ([31$ and the facts that 2||/||ooe' = e and 



(31) 



we obtain that 



P'^ I max sup 



I E {/(C^/CD-/(C?,,,,^,)_,/CJ,,,^,^,)J} >e)<C 



lan^e-^"'^ (32) 



r-t,<n 



for some constants C12, C13 > 0. 

Let us now consider the first term in the r.h.s. of (|27p . which corresponds to the evolution due 
to mutation. Recall that B denotes the generator of the mutation process, and that mutations occur 
independently along distinct lines. Recall also that we write Bf for the action of B on the second 
argument of /. We have 



j-lj<n 



+ - E 



j:ij<n 



Bf{aM)ds 



1 ff'^'^ij 



ds, 



j:ij<n 



where (Aj(f))(g[^. q.^^j is a zero-mean martingale (as the sum of finitely many zero-mean martingales). 
Recall the event An defined in (|25|) . For any n large enough so that 29n~'^\\Bf\\oo < e, we have 



P I max sup 



^ E i^^^(tATij)-^^{tATij)-) /(Ca,'^a»)} 



j-Ij<n 



> 2e 



< P-(^„) + VP^f sup |A,(t)| + ^n-2||i3/|U > 2e; aA 
<P-(A„) + EP"f sup |A,(t)|>e;^J^y 

Hence, using Lemma 3.1(6) in |BB+09| together with the same argument as in (j26p to control the 
number of individuals in Sf and with label at most n at any given time, we obtain that 



P'^ I max sup 



\ E >25) <(2(TT + r)n2 + l) 



-Ci4n 



j-Ij<n 
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for some C14 > 0. Together with (I27p and (j32p . this gives us that 

P'^f max sup |(M«,/) - (M^,,/)| > 3e) < Ci^n^-^^^^ (33) 
for some C15, Cie > 0. 

Now, it is easy to see that if maxi<j„ Hi < e and maxj<i„ sup„^<j<^^^j |(M;', /) - {M^^, f)\ < 3e, 
necessarily dj > Oj+i for every i, which imphes that 

max sup |(Mr,/)-(M^,/)| <6e 
on this event. Using (j26p and (j33p . we finahy obtain that 

P'^f sup |(M-,/) - (M-,/)! > lie) < Cun^e-^i^'^ =: 5„, 

\0<t<T / 

and Lemma 14.51 is proved. Then, the Borel-Cantehi lemma gives us the a.s. convergence of (M",/) 
towards (M°°,/) as n — t- 00, uniformly over compact time intervals. Since by assumption the set 
V^B) n Cc(T X K) is separable and is dense in Cc(T x K), the a.s. convergence of M" towards M°° 
follows from Lemma 12.21 as well as the fact that M°° has cadlag paths with probability 1. 

Finally, let us now prove that M°° has the same law as the quenched spatial A-Fleming-Viot process 
of Theorem 1. As we have already mentioned at the beginning of Section [4.21 with the notation of the 
proof of Theorem 1' we have 

n 

Letting n tend to infinity and using the a.s. convergence of n^^M^ towards the random measure 9Ko,t, 
we can conclude that has law Qq^^iti, .). 

Let us now fix < s < t and condition on = m € Mx- By construction, the types of 

all the individuals alive at time t are determined by propagating the types of the individuals living 
at time s along the genealogical trees created by the look-down dynamics. But Lemma 14.21 and the 
Markov property of the genealogical processes ensure that the ancestries between times s and t of any 
finite sample of individuals living at time t have the same law as the trees {Ah)o<h<t-s under the 
probability measure P"^'* (recall that A is defined in the paragraph around ([2]), and P'^'* is introduced 
just after Lemma l2.ip . Consequently, the proof of Theorem 1' shows that has law Qg^{m,d'm'). 
We can therefore conclude that M°° is a Markov process with transition semigroup {Q'^ t)o<s<t, and 
so Theorem 1 guarantees that M°° has the same law as the quenched SLFV {Mt)t>o- The proof of 
Theorem 2 is now complete. □ 

Let us end this section with the proof of Lemma 12.61 

Proof of Lemma 12.61 Again, let us consider / G Cc(T x K) and show that {{Mjr^, f))t>o has paths of 
finite variation with P'^-probability 1. By taking a countable basis /i, /2, . . . of Cc(T x K), we shall then 
conclude that this property holds for all /j's simultaneously with probability 1, and so by Lemma 12.21 
that a.s. M°° has paths of finite variation in Dj^^ [0, 00). 

Let us thus fix T > 0. Since there is no mutation, M°° evolves only at the times of a reproduction 
event. By construction, during the event (t,z,r,u), at each site x within M{z,r) a fraction u of the 
population is replaced by individuals that are all of some type k. That is (in the notation of ([7])), 
pf{x) = (1 —u)p'^{x) +u6k- Together with the fact that the spatial marginal of each Mj°° is Lebesgue 
measure, for any finite number of times to, . . . ,ti such that = to < ti < . . . < ti = T , we have 
l-i 

j=0 {t,z,r,u)euj:0<t<T 

<2||/|U E uYol{Sfnn{z,r)), 

(t,z,r,u)eiu:0<t<T 
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where here again Sf stands for the compact support of / with respect to the spatial coordinate. But by 
(j23p and (j24p . the sum in the r.h.s. above is finite for P-a.e. lo, and is independent of the subdivision 
(to, ■ ■ ■ ,ti) chosen. Since the total variation of ((Mj°°, /))o<t<T is given by 



sup sup 5;|(M,-^,/)-(M,7,/)| 

«6N (to,...A)^ 



the desired result follows. □ 
4.3 Generalizations 

The construction carried out in the last paragraph is robust in the sense that it can accommodate 
different variants of the spatial A-FIeming-Viot process such as that introduced in |BKE10| . In the 
Gaussian model of ^BKElOj . the environment is given by a Poisson point process of epochs and centres 
with intensity measure cdt® dx, where c > 0. During an event centered on z € T, the fraction u{z, x) 
of the population killed at site x is given by 



where uq G (0, 1] is the maximal killing intensity and 0^ > is a fixed parameter. Then, the location 
of the parent is sampled according to the kernel 



where a > (note that if a > 1, the parent is chosen from a wider area than that over which the 
impact of the event is non-negligible). A type is then drawn from the local type distribution at this 
site. Finally, the population removed during the first step is replaced by offspring of the chosen parent 
so that the local mass of individuals at every site remains constant. This process can be phrased as an 
A^A-valued evolution, and the corresponding genealogical process can be described explicitly. 

In order to adapt the look-down construction of Section 14.11 let us emphasize its two main ingre- 
dients: 

An ancestry given by a finite-rate jumps process: A lineage at location x is affected by an event 
whenever it belongs to the fraction u{zi,x) replaced during the event. Hence, with (|34p . the rate 
at which a lineage jumps is equal to 



This fact enables us to introduce a well-defined procedure to make individual types evolve in 



Reversibility of Lebesgue measure: With the rule (|35p . the transition kernel Ki corresponding to 
the i-th. event, centered on Zi, is given by 



Since conditionally on Zj the new location of C is independent of its previous position, it is not 
difficult to see that 

dxKi{x,dy) = dy Ki{y,dx). 
Hence, Lebesgue measure is again reversible for the evolution of C- 




(34) 




(35) 




time. 



Ki{x,dy) = (1 - u{zi,x))6x + u{zi,x) v{zi,y) dy. 
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Next, when a = 1 we can again choose the parent to be the affected individual with lowest level. When 
a ^ 1, the kernel with which the parent is drawn differs from that with which individuals are affected, 
and so we have to change the way the parent is sampled. This time, we thin the Poisson point process of 
individuals at the time of the event by keeping each point with probability C|.„)/t'(0, 0). 

We then take the individual of the thinned point process with lowest level as the parent, and resample 
its location using the kernel v{zi,y) dy. In this way, some of the affected individuals may now have a 
lower level than the parent. However, the resampling of the location of the parent ensures that the 
position of the ancestor just after a given merger of V has density v{zi, ■ ), which is the essential point 
in the equality in distribution of V and A (cf. Lemma l4.2p . 

Using the two properties emphasized above and the ideas of Sections 14.11 and 14.21 one can then 
construct a particle representation for the stochastic flow obtained in the Gaussian model. 

Another possible generalization is to allow multiple parents. Extending the ball model from Sec- 
tion 14.11 we may e.g. fix a distribution 7 on N with compact support and associate an independent 
realization Ni of 7 to each event of lo. During the i-th event, Ni parents are picked independently 
and uniformly over the area of the event. If the corresponding types are ki, . . . , k^. (not necessarily 
distinct), the new value of the SLFV is given by 




In words, within the area of the event we keep a fraction 1 — Uj of the population as it was just before 
the event, and replace a fraction Ui by offspring of the Ni reproducing types in equal proportions. 
Thus, the associated genealogies can have multiple and simultaneous mergers whenever 7 puts some 
mass onto {2, 3, . . .}. 

In the corresponding look-down construction, the only difference with that of Section 14.21 is that 
we use the Ni affected individuals with lowest levels as the parents during event i, and decide that 
an affected individual looks down onto one of these Ni levels with equal probability, independently of 
each other. Since the two conditions given above are fulfilled (with the kernel Ki given by (|20p ). we 
can construct a look-down coupling between the quenched SLFV and its genealogies. 

5 Coming down from infinity 

Recall that a coalescent is said to come down from infinity if, starting from countably many lineages, 
there exists a time in the past at which the number of ancestors is finite. For non-spatial exchangeable 
coalescents, it is known that whenever the quantity corresponding to the impact u here is always less 
than 1, then either the coalescent comes down from infinity instantaneously with probability 1, or the 
number of ancestors remains infinite for all times a.s. See e.g. Proposition 23 in |Pit99| for a statement 
of this result for coalescents with multiple mergers, and Lemma 31 in [SchOOj for the more general 
case of coalescents with simultaneous and multiple mergers. Furthermore, a precise criterion for GDI 
is obtained in |SchOOb| for A-coalescents, and some conditions are given in Section 5.5 of |SchOO| as 
concerns the H-coalescents. 

In the context of spatial coalescents, only a few results exist concerning the question of coming 
down from infinity. Indeed, the geographical movement of lineages may separate them before they find 
a common ancestor, or may bring back together some lineages which were too far away from each other 
to coalesce. Hence, understanding the form of the resulting genealogical process requires a fine analysis 
of the interplay between these two mechanisms. In [LS06] , Limic and Sturm consider a population 
spread over the vertices of a finite graph. They assume that the lineages migrate independently of 
each other between the sites and can coalesce only when they belong to the same subpopulation. They 
show that the corresponding spatial coalescent comes down from infinity if and only if its non-spatial 
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counterpart does. The case of discrete but infinite graphs is then explored in |ABL10| . where it is 
shown that because the timescales of migration and coalescence are precisely the same, an arbitrarily 
large number of lineages can escape from the others without coalescing. As a consequence, the spatial 
coalescents they consider never come down from infinity. Finally, a version of the SLFV in which the 
geographical space is a self-similar hierarchical structure (such as a Cantor set, or an infinite m-ary 
tree) is introduced in |Frell) . The forwards-in-time evolution is formulated as a stochastic flow in 
the same spirit as Bertoin and Legall's construction of the A-Fleming-Viot superprocess in |BLG03| . 
Freeman then details the five different possible behaviours of the corresponding ancestral process, 
making explicit the importance of the geographical structure in the form of the genealogies. 

Coming back to the SLFV, imagine we sample countably many individuals in the ball B := B(0, 1) C 
T and trace back their genealogical relations. We want to show that, under the condition 



CXD 



z.,({l})/i(dr) =0, (36) 







in a.e. environment the number of distinct ancestors remains infinite for all times in the past. 



Remark 5.1. The condition l[36\) simply says that the impact parameter is always less than 1. If it 
does not hold, an event overlapping all the lineages but for a finite number ( say, n) of them may have 
impact 1, in which case the number of ancestors would come down to n + 1 < oo. 

In fact, as one might expect from the comparison of the finiteness condition ([T]) with the condition 

Rate at which a sing le lineage is affected = / ^ < oo 

Jo u 

for the presence of dust in a nonspatial A-coalescent (see Theorem 8 in pit99j ). we shall show that at 
time t > in the past infinitely many lineages have not yet been affected by an event. The difficulty 
here comes from the spatial correlations between the rates at which close-by lineages are affected by 
reproduction events. Despite this effect of space, the proof of the following proposition relies on the 
Poissonian structure of the events. Recall the condition ([T| of existence of the spatial A-Fleming-Viot 
process on T with parameters ^ and {ur, r > 0}. 

Proposition 5.2. Assume that /i36\) holds. Let M be a Poisson point process on B x [0, oo) with 
intensity measure dx\^ (8) dl, and let us use {({«}, Xj) : {xi,£i) € M} as the initial value of the spatial 
A-coalescent A defined in Section \2.1\ Then for almost every environment oj, at any time t > the 
set of the ancestral locations in At contains a Poisson point process on B with infinite intensity. In 
particular, the spatial A-coalescent A never comes down from infinity. 

Proof of Proposition 15.21 For every x B, let us write a{x,t) for the set of events of u between 
(backward) times and t overlapping x. We have 

P(P^'°(C does not jump until time t) > O) = P( JJ (1 - Ui) > ) = 1. (37) 

^ i€a{x,t) ^ 



explt / Vr,x [(1 - - 1] iJr{du)fi{dr) \ , (38) 



Indeed, using the exponential formula for the Poisson point process of events, we can write for every 

( fOO f\ 

]E EiSa(a;,t)log(l-«i) 

'0 ^0 

where Vr,x stands for the volume of the ball i?(x, r) in T. But for every u G [0, 1), 1 — (\ — uf decreases 
to as decreases to 0. Furthermore, using this property with ^ = 1 we obtain that for every € (0, 1), 

0< / / Vr.x\^ — — u)^\vr{du)^{dr) < / / Vr^xUVr{du)[l{dr) < (X) 

Jo Jo Jo Jo ' 
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by the condition ([T]) imposed on and {ly^-, r > 0}. We can therefore use dominated convergence and 
the fact that the set {u = 1} is never charged to conclude that 



lim 



oo />1 



Vr,x [(1 - uf - l] Vr{du)n{dr) = 0. 



Coming back to (j38p and letting 6 tend to 0, we obtain that lP'[X^i£cj(a; t) log(l 
thus (|37p is proved. 

Using now Fubini's theorem, we can write 



Ui] 



-ool = and 



= E 



\o\{B) Jb "^{n,6«{.,t)(i-«.)>o} "^"^ 



and since the quantity within the expectation on the right-hand side belongs to [0, 1], it therefore needs 
to be 1 P-a.s. 

In particular, for a.e. u, we can find £{uj) € (0, 1) such that 



Let us call Se^ui the set {x £ B : Y\i^ct{x t)(^ ~ ^ ^(^)}- What we have just shown is that for 
almost every ^s^lj has positive Lebesgue measure. Hence it contains infinitely many points of A/^, 
and each of them has probability at least e(w) not to have been affected yet, independently of each 
other. Together with initial Poissonian structure of J\f, this completes the proof of Proposition 15.21 □ 
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